% A = [1,2,4;
%     0,1,1;
%     0,0,9];
% b = det(A);
% c = inv(A);
% b
% c
% fun = inline('exp(x).*sin(x).^2');
% % quadl(fun,0,1);
% quadl(inline('exp(x).*sin(x).^2'),0,1)
% taylortool
% det([1 2;3 4]+10-2i) 1：一个矩阵，该矩阵每个位置加上一个数字，该数字为虚数
% [1 2; 3 4].*[0.1 0.2; 0.3 0.4]  两个矩阵，对应的位置相乘
% [1 2; 3 4].\[20 10;9 2]         两个矩阵，对应的位置相除，后者除以前者
% [1 2; 3 4].^2                     矩阵的每个位置开平方
% exp([1 2; 3 4])                     矩阵每个位置的值以e为底求方
% log([1 10 100])                     求矩阵对应位置的值是e的几次方
% prod([1 2;3 4])                       求矩阵每一列元素相乘的积
% [a,b]=min([10 20;30 4])                 求矩阵每行中最小的元素，第一个参数为最小的值，第二的参数为对应的下标
% abs([1 2;3 4]-pi)                   矩阵上的值减去 pi 让后新矩阵求绝对值
% [10 20;30 40]>=[40,30;20 40]       矩阵对应的值比大小，小的话为0否则为1
% find([10,20;30,40]>=[40,30;20,10])      比较对应的位置的大小，把位置列成线性，返回为真的下标
% [a,b]=find([10 20;30 40]>=[40,30;20 10]) %(提示：a为行号，b为列号)  比较对应的位置的大小，返回为真的下标
% all([1,2;3,4]>1)          判断每个子元素是否都为真
% 
 % any([1 2;3 4]>1)         判断每个子元素是否有一个为真
% 
% linspace(3,4,5)          从3开始到4结束,平均分配五个值
% 
% A=[1 2;3 4];A(:,2)      从头开始到尾结束,步长为2遍历矩阵



% % clear; a=1,b=num2str(a),c=a>0, a==b, a==c, b==c
% clear; fun='abs(x)',x=-2
% b = eval(fun)
% double(fun)

% r = 2;
% p = 0.5;
% n = 12;
% 
% T = log(r)/(n * log(1 + 0.01*p));
% fprintf('%d',T);

% z=magic(10);
% 
% ret2 = sum(z);
% ret3 = sum(diag(z));
% z(:,2) = z(:,2)/(3^(1/2));
% z(8,:) = z(8,:) + z(3,:);
% b = Untitled;
x= [81, 70, 65,51, 76, 66, 90, 87, 61, 77];

pos = length(x);
numsum = 0;
while(pos)
    numsum = numsum + x(pos);
    pos = pos-1;
end
ret1 = numsum / length(x);
fprintf("%d",ret1);







